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Statics of polymer droplets on deformable surfaces 
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Friedrich-Hund-Platz 1, 37077 Gottingen, Germany 

The equilibrium properties of polymer droplets on a soft deformable surface are investigated by molecular 
dynamics simulations of a bead-spring model. The surface consists of a polymer brush with irreversibly end- 
tethered linear homopolymer chains onto a flat solid substrate. We tune the softness of the surface by varying 
the grafting density. Droplets are comprised of bead-spring polymers of various chain lengths. First, both 
systems, brush and polymer liquid, are studied independently in order to determine their static and dynamic 
properties. In particular, using a numerical implementation of an AFM experiment, we measure the shear 
modulus of the brush surface and compare the results to theoretical predictions. Then, we study the wetting 
behavior of polymer droplets with different contact angles and on substrates that differ in softness. Density 
profiles reveal, under certain conditions, the formation of a wetting ridge beneath the three-phase contact line. 
Cap-shaped droplets and cylindrical droplets are also compared to estimate the effect of the line tension with 
respect to the droplet size. Finally, the results of the simulations are compared to a phcnomcnological free- 
energy calculation that accounts for the surface tensions and the compliance of the soft substrate. Depending 
on the surface/drop compatibility, surface softness and drop size, a transition between two regimes is observed: 
from one where the drop surface energy balances the adhesion with the surface, which is the classical Young- 
Dupre wetting regime, to another one where a coupling occurs between adhesion, droplet and surface elastic 
energies. 



I. INTRODUCTION 

The molecular structure and dynamics at the bound- 
ary between a liquid and a confining surface dictate the 
interactions of a liquid with the surface and the friction 
of a fluid that flows past iti"— . They determine the equi- 
librium shape of drops and the microscopic dissipation 
mechanisms give rise to slippage, control the (de)-wetting 
kinetics of thin coating layers^"— and set the time scale, 
on which droplets spread on a wettable surface. Since 
early works of YoungH and Navier— , the effect of molec- 
ular interaction and friction at the surface are described 
within a continuum theory via boundary conditions, e.g. 
Young's equation and Navier's slip condition, featuring 
the surface and interfacial tensions, the hydrodynamic 
interface position, and the slip length as coarse-grained 
parameters that transfer the microscopic information to 
the continuum description. 

Which coarse-grained parameters describe the micro- 
scopic behavior depends on the specific system. Often, 
idealized hard planar solid surfaces are considered. Such 
ideal surfaces, however, are rarely encountered in prac- 
tice, and the softness of a surface is a natural departure 
from the ideality of surfaces. These deformable surfaces 
are expected to differ in their static wetting behavior and 
the softness also influences the flow of liquids past it. 

For instance, when a droplet is deposited on a sur- 
face, and the spreading process has been completed, an 
equilibrium contact angle Oy is reached. Classically, this 
spreading process should essentially depend on the capil- 
lary properties of the liquid and the surface, and also on 
the liquid viscosity. The equilibrium contact angle for a 
non-deformable surface, is given by the parallel balance 
of surface-tension forces at the contact lineH. For soft 
elastomeric surfaces, however, a local deformation will 



occur near the triple line, which is dictated a balance be- 
tween the component of the surface tension of the liquid 
perpendicular to the surface and the ability of the soft 
surface to deform. This interplay gives rise to a lifting-up 
of the contact line, which has been shown to largely con- 
tribute to the wetting or dewetting on soft elastomeric 
surfaces^ ^^^'^"" — , as well as to the visco-elastic braking 
of a driven droplet rolling over a deformable surface^. 

Molecular dynamics (MD) simulations provide a well 
suited framework to investigate the physical quantities 
pertinent to the wetting of liquids on soft deformable 
solid surfaces. In the following, we study a polymer 
brush as a prototypical example of a soft deformable sur- 
face, where flexible polymers are tethered with one end 
to a flat, impenetrable and hard substrate. The ability 
of the soft surface to deform is controlled by the grafting 
density of tethered chains. The softness of the surface 
is, in a sense, universal because it stems from the bal- 
ance between the loss of conformational entropy as the 
tethered macromolecules stretch away from the grafting 
substrate and their excluded volume interactions. Both, 
in our simulations as well as in experiments, it can be 
tuned by altering the degree of polymerization of the 
brush molecules or their grafting density without chang- 
ing the non-bonded interactions. This soft surface is in 
contact with a polymer liquid. Varying the incompati- 
bility between the brush and the polymer liquid we can 
independently tune the softness and the contact angle 
that a liquid drops forms on the brush. 

Our work is organized as follows: In sec. |II] we de- 
scribe the model and provide details of the simulations. 
In the Sec. Ill A[ the equilibrium properties of droplets 
properties are studied, and the deformability of the soft 
surface is characterized in the Sec. IIIBI The observa- 
tion of the equilibrium wetting of cylindrical droplets in 
contact with different deformable surfaces are discussed 



in the Sec. IIII Al Then, the effects of the hne tension 
are investigated in the Sec. IIII Bl in which spherical, cap- 
shaped and cyhndrical droplets in contact with the same 
soft surface are considered. Finally, we summarize our 
conclusions and give a brief outlook. 



II. MODEL AND SIMULATION ASPECTS 

MD simulations are performed using a well established 
polymer coarse-grained modeled. Each polymer is repre- 
sented by a chain of N beads of mass m = 1 connected 
by springs to form a linear, flexible chain. Beads, that 
are not neighbors along a chain molecule, interact with a 
Lennard- Jones (LJ) potential: 
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for r < Tcut = 2.5(7, while for r > Tc, the potential is 
cut-off and shifted such that it is continuous at rcut ■ The 
potential comprises a steep repulsion at short distances, 
which mimics the excluded volume of the beads, and a 
longer-ranged attraction, which makes the polymer con- 
dense into a dense liquid that coexists with its vapor of 
vanishingly low density (poor solvent condition). Indices 
(a, 13) stand for the different types of pairs, dd, bb, and bd 
for droplet (d) or brush (6). All energies are measured in 
units of e. Adjacent beads along the chains are coupled 
through an anharmonic finite extensible nonlinear elastic 
potential (FENE): 
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Table I. Bulk properties of the polymer liquid at fcsT = 1.2e 
and pi ^0.79a'^. 



where f,;j = (r^ — Yj)/rij and v.^ 
the weight functions: 
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with Tc identical to the one used in Eq.([T]). The random 
force is given by: 



Ff = ^^..«(r,,)%f,, , 
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where 9ij is a random variable with zero mean, 
unit variance, second moment {9ij(t)9ki{t')) = 
(SijSji + 6iidjk)6{t - t'), and % = %. The weight 
functions Lj^{rij) satisfy the fluctuation-dissipation 

theorem, [w^l = uj^ . Friction F and noise strength 
£_ define the temperature via £,^ = 2kBTT. We choose 
r = 0.5r~^ in all our simulations. 
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where model parameters are identical to those given in 
Refi^, namely k = SOe/a^ and i?o = l.Scr, chosen in 
such that unphysical bond crossings and chain breaking 
are eliminated. The bonded interaction, Eq. ([2]), is ap- 
plied in conjunction with a purely repulsive LJ potential, 
i.e. Eq. ([1} with rcut = 2^/^cr. All quantities are ex- 
pressed in terms of molecular diame ter a ~ 1, LJ energy 
e = 1, and characteristic time r = -y/mcr^/e. 

The classical equation of motion are integrated via 
the velocity- Verlet algorithm with a time step of Ai = 
5 X 1Q~^T. Temperature is kept constant at ksT = 1.2e 
using a dissipative particle dynamics (DPD) thermostat, 
that also accounts for hydrodynamic interaction due to 
the local conservation of momentu m^^i^^ . The thermo- 
stat adds to the total conservative force, that arise from 
([1]) and ([2]), a dissipative force Fj on each monomer i, 
and a random force, Yf. Both forces are applied in a 
pairwise manner, such that the sum of thermostatting 
forces acting on a particle pair vanishes. Let F be the 
friction constant, the dissipative force is given by: 
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A. Properties of the polymer liquid of the droplet 

Liquid droplets consist of M = 450, • • • , 15328 poly- 
mer chains of polymerization degree of A^ = 10, 100 and 
200 beads. Only the two last degrees of polymerization 
may involve entanglement effects^, as the entanglement 
length is A^o ~ 64 for this model. The interaction param- 
eter in Eq. (jlj was set to Cdd = l.Oe for the polymer liquid 
of the drop. For the temperature ksT ~ 1.2e, the coexis- 
tence density of the fluid with its vapor is p\ ~ O.TQc"^, 
while the vapor density can be negligibly small^*^. 

Bulk properties of droplets are summarized in Tab. HI 
The Rouse relaxation time Tp is related to the self- 
diffusion coefficient D = R'^/ (Sir^Tp), where Re = 
(i?c)^^^ is the average end-to-end polymer distance. For 
the smallest chain length, A^ — 10, this corresponds to 
very small values of the invariant degree of polymeriza- 
tion, Af = {pxRllNf w 14. 

The shear viscosity, 77^, has been computed using the 
reverse non-equilibrium molecular dynamics method^i^, 
which provides an efficient means compared to other 
(non)-equilibrium molecular dynamics methods^. To this 
end, one divides the system into slabs and exchanges 
momenta between the two slabs that are separated by 
half the system size. This creates a momentum flux 




Figure 1. (Color online) Sketch of the simulation set-up at 
ksT = 1.2e. Brush segments are colored in dark green, 
grafted ends are depicted in light green, while liquid droplet is 
colored in blue. Periodic boundary conditions are used in both 
X and y direction. The figure corresponds to a snapshot of the 
simulation box for a brush of grafting density pgRc = 9.5, and 
chain length A*' = 40, while the droplet contains M = 1916 
chains of length TV = 10. 



between slabs and results in a linear velocity profile 
v\i{x). By varying the number of momentum swaps 
Ns = 3,15,60,120,300, or 1200 per integration step, 
one tunes the effective shear rate imparted onto the sys- 
tem. Simulation runs, which lasted between SOOOr and 
24000t, were performed in a simulation cell containing 
19160 beads. Data were collected at half these times, 
when a steady state has been reached. The values ob- 
tained using this method are then extrapolated to van- 
ishing momentum flux. Results are given in Tab. U and 
are found to be in agreement with ones, for example, from 
equilibrium molecular dynamics simulations^i^^. 

Additionally, we computed the surface tension 7d, 
using a slab geometry^^i^i. The anisotropy of 

the pressure tensor yields the estimate 2A'yd = 
V {{Pzz) — 0.5 {{Pxx) + {Pyy))}i where the cross section 
area is ^ = L^; x L^, and the factor of 2 arises from the 
number of Hquid/vapor interfaces in the simulation box 
with periodic boundary conditions. 

Except for Sec. IIIIBl cylindrical droplets are consid- 
ered with a symmetry axis in the x direction, along which 
periodic boundary conditions are applied. A snapshot of 
such a system is shown in Fig.[TJ In the following, brushes 
with an aspect ratio A^ = Ly/L^ = 5 will be used, also 
with periodic boundary conditions in y direction, while 
impenetrable, flat surfaces limit the simulation box in 
z direction. For large drops and small values of L^, 
the cylindrical state is then an absolute minimum of the 
surface free energy, i.e. such a geometry stabilizes the 
cylindrical droplet against cap-shaped droplets, which 
will become energetically favorable for small droplets or 
larger system sizes, L^- Cylindrical droplets are charac- 



terized by straight contact lines. Thus, we avoid strong 
finite-size effects which arise for spherical drops due to 
dependence of the length of the contact line on the 
droplet size and the concomitant contribution of the line 
tension^i^i^i^i^. 



B. Characteristics of the soft, polymer brush 

The deformable surface is modeled by an end-grafted 
brush of polymers comprised of A^ = 40 beads. The inter- 
action parameter in Eq. ([Ij is identical to the one used for 
the liquid droplet, namely ebb = l.Oe. The grafting sub- 
strate is hard and impenetrable and the van-der-Waals 
interactions between the solid and the polymer brush and 
liquid are catered for by an integrated Lennard- Jones po- 
tential: 
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-const for z < 2.5cr 

(6) 

^waii = for z > 2.5(7, and the constant in Eq. dH) is 
chosen such that C/]^j'" is continuous. We set Cwaii = 0.0 le 
and a = 1.0. Since e^aii 'C fcsT, the wall substrate is 
effectively repulsive. 

The first anchor segment of each polymer of the 
brush is irreversibly grafted and immobile. The an- 
chor beads are regularly placed on a square grid in the 
xy-plane at a distance zq = 0.8a above the substrate. 
The bond, Tq — ri, between the immobile, anchor and 
the second, mobile bead of the polymer brush is rep- 
resented by an anharmonic spring of the form U^ ~ 



Ea (ri - ro) • A^ - (ri - ro)' 



27 



with 



1.5e and 



A = 0.8(7. This anharmonic potential is softer than the 
FENE potential that is used for the other bonds and thus 
copes better with the strong forces on the bonds at the 
grafting substrate. The grafting density, pg is defined by 
the number of polymer chains per unit area. We measure 
it in units of the end-to-end distance Re = 7.74(t of melt 
chains of length A = 40. A value of p* = pgi?^ - 0{1) 
marks the crossover from mushrooms to a proper polymer 
brush. Systems with reduced densities of 6.5 < p* < 57 A 
were investigated. This wide range of grafting densities 
represents different physical regimes for the deformable 
surface: A low grafting densities the surface is very 
soft; at high densities, the brush polymers are strongly 
stretched and the surface is significantly less deformable. 
First, we study the height, H, of the brush as a func- 
tion of the grafting density pg. The average height, H, 
is then defined as the first moment of the density distri- 
bution H = J z^{z)dz/ J ^{z)dz. We confirm that this 
quantity scales like H ~ vNp^, where v is the effective 
volume of a coarse-grained segment. The exponent, x, 
depends on the solvent quality. Since the temperature 
ksT = 1.2e corresponds to a polymer brush in a poor 
solvent, the brush is nearly incompressible and i; ^ 1/pi- 
Thus, the average height of the brush scales linearly with 



grafting density, x = 1, in agreement with previous sim- 
ulations^. 

Under bad solvent conditions, the top surface of a 
brush resembles the interface between a dense poly- 
mer liquid and its coexisting vapor of vanishingly small 
density. Conceptually, one has to distinguish between 
the "intrinsic" width of the brush surface and long- 
wavelength fluctuations of the local height of the brush, 
which are the analog of capillary waves on a free liquid- 
vapor interface. The long-wavelength fluctuations of the 
local height of the brush can be described in terms of an 
effective interface Hamiltonian^: 



■Hcap[^int(a;,2/)] 



d^r. 



[7b [VZint]^ + 

Kh {Zint - Zintf) , (7) 



where r|| = {x, y) denotes the two lateral coordinates 
parallel to the grafting surface of area, A — L^Ly. The 
variable Zintix,y) is the instantaneous local position of 
the brush- vapor interface, and Zmt denotes the average 
height H. 7^ denotes the surface tension of the brush and 
Kfi characterizes the coupling between the average brush 
height and the grafting substrate. Since the brush is a 
slightly compressible liquid, fluctuations of the average 
brush height occur and their free-energy cost is propor- 
tional the the area of the grafting substrate and scale 
quadratically with the deviation from H. It is useful to 
consider the Fourier components of the local interface 
position. 
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with wave vectors q^ = 2t: {kx/ L^yky/ Ly) and kx,ky G 
N. With this coarse-grained description in Eq. (O, 
the equipartition theorem connects the spectrum of the 
brush- vapor interface fluctuations to the tension and de- 
formability of the surface: 



(c^(q)) = (ciiq)) 



2kBT 



A (7bq^ + Kb) 



(9) 



In order to evaluate Eq. ([9]), the brush surface is 
divided into a grid of columns with lateral extension 
A^ ~ Ay ~ 1(7. In each column, (x, y) the brush- 
vapor interface position Zint{x,y) is defined as a Gibbs 
dividing interface^i^ between the vapor and the brush 
melt. Data in a distance of 4cr away from the local in- 
terface position have been analyzed to minimize the ef- 
fect of bulk-like density fluctuations of the compressible 
brush. About 10"^ configurations with an aspect ratio 
Ar = 1, separated by 5r, were analyzed. The roughness 
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Figure 2. (Color online) Surface tension 7b of the A'' = 40 
beads per chain deformable surface as a function of the re- 
duced grafting density Pg — p^Rc, extracted from the fluctu- 
ations spectrum, Eq. ((9]). The horizontal dashed line is the 
liquid-vapor surface tension of a polymer melt with identical 
chain length. Inset: Elastic constant Kb of the brush surface 
as a function of the reduced grafting density. The data are 
compatible with a linear increase with the reduced grafting 
density as indicated by the dashed line. 



average is performed over all grid columns and time, was 
found to decrease (not shown) with an increase of p*, as 
expected when the brush stretches. 

The Fig. [2] shows the values of the normalized brush 
surface tension j\ya^ /kBT (main panel) and brush's elas- 
tic constant k\,(t^ /ksT (inset), extracted from the small q 
behavior of the y-averaged Fourier components Cc{qx) = 
Cc{qx, % = 0)1 as a function of the reduced grafting den- 
sity. For small grafting densities, the tensions 7b in- 
creases and saturates at larger p*. Indeed, at interme- 
diate grafting densities, the structure of the brush in a 
bad solvent is independent^^ from pg and its surface re- 
sembles a liquid-vapor interface. Thus, we expect that 
the surface tension, 7b, of the effective interface Hamil- 
tonian is comparable to the surface tension of the free 
polymer liquid in coexistence with its vapor. For the pa- 
rameters used in our simulation. fc^T = 1.2e and A^ = 40, 
the liquid-vapor interface tension is 7dcr^/e = 0.67 ±0.02 
and this value is indicated by the horizontal line in Fig. [21 
Indeed, our simulation data are compatible with 7b — > 7d 
in the limit of large grafting densities. 

On the other hand, the clastic constant Kb, which is 
related to a uniform height variation, appears to de- 
pendent more strongly on the grafting density and the 
simulation results are compatible with a linear increase 
with the brush grafting density. If one assumed that 
each chain acts like an entropic spring with an elas- 
tic constant fc, which is a reasonable assumption for a 
brush in contact with a melt of identical polymers^, that 
Kb — BfcsTpJ/i?^. For our system, we observe the same 




Figure 3. (Color online) Sketch of the simulation analog of an 
AFM-indentation experiment and illustration of the quanti- 
ties related to the measurement of the effect shear modulus, 
Gb, of the brush. In the direction of the applied force, the 
brush responds like a liquid, while due to the constraints im- 
posed by the grafting, the brush's response shows a resistance 
to the resulting lateral force. 



linear scaling with the grafting density but the pre-factor 
is about a factor 5 larger. This difference presumably 
stems from the fact that our system is much less com- 
pressible than the brush-melt system studied in Ref4^. 
In the brush-melt system the average height of the brush 
can fluctuate without a change of the segment density 
inside the brush - free polymers of the melt in contact 
with the brush can enter or leave the brush to compen- 
sate for the variation of the density of brush segments. 
This mechanism is not available for the polymer brush in 
contact with vacuum, where a fluctuation of the average 
brush height entails a change of the polymer density and 
a concomitant free-energy penalty due to the equation- 
of-state of the polymer liquid. 

Alternatively, the deformability of the surface can be 
characterized by its shear modulus, Gb- Using computer 
simulations, this quantity is often difficult to access, while 
experimentally this quantity can be efficiently measured 
by indenting an AFM tip into the soft surfac o^^'^^ . Pre- 
vious theoretical work by Frcdrickson et alA^ predicted 
that a molten polymer brush has an effective shear mod- 
ulus, Grmb, that causes the indentation behavior of a soft 
brush surface to deviate from that of a liquid. The shear 
modulus can be estimated by applying a lateral generic 
force, F, to the top of the brush surface, which "shears" 
the chains by displacing each chain end at the surface 
by a distance x, as depicted in Fig. |31 Let Atip be the 
lateral contact area of the indenter and ctip the depth 
of the indentation. Following the qualitative considera- 
tions of Ref4^, one estimates that the free end of a typ- 
ical brush molecule, which is located a distance y/Aup 
away from the center of the indenter, is laterally dis- 
placed by an amount x ~ i/pg^tip x ^tip/{y^H), where 
the first factor accounts for the distance of the chain 
from the center of the indenter in units of the distance 
of the grafting points. The second factor is the lateral 



displacement of a chain at the center of the indented re- 
gion. Such a lateral displacement gives rise to a lateral 
force, fx ~ ksTx/Rg, on each chain and the total lateral 
stress generated by all chains underneath the indenter is 
<^xz ^ Pgfx ~ PghsTx/R^. Since the trajectories of the 
chains are approximately linear, the strain Uxz is of the 
order x/ H . Conceiving the brush as an incompressible, 
isotropic, elastic film^, one can define an effective shear 
modulus of the brush via the relation axz = GhUxz, i-e. 
Gb ~ kBTpgH/R^. Using H ~ vNag and including the 
numerical pre-factor of the more detailed calculation in 
the framework of the Alexander-brush model^, one ar- 
rives at: 
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kBTvN ,2 

i?3 i?3 ^S 



(10) 



Only recently, Fujii et a/.— have shown, by measuring 
the power spectral density of monodisperse polystyrene 
brushes submitted to a tapping-mode AFM tip, that 
such theoretical predictions can describe the experimen- 
tal data. 

As illustrated in Fig. [3l the AFM-like method is nu- 
merically mimicked by forcing a spherical indenter of ra- 
dius, R, to compress the deformable surface. The lat- 
eral coordinates, (xtip,2/tip), of the indenter are fixed. 
The spherical indenter exerts a radial force of magnitude 
F{r) ~ ~K{r — R)^ on each segment, where the force 
constant is K = 100 and r stands for the distance of 
the segment to the center of the indenter. The force is 
repulsive for r < R, and F{r) = for r > R. In the fol- 
lowing, we use large tips, whose radii are greater than the 
brush height. This set-up corresponds to a compressive 
mode of indentation. Two types of tips are considered: 
(i) a tip with a fixed radius R = Raxed and a center 
of mass Tj^ (t) that weakly varies in time as to describe 
the approach of the tip towards the brush surface, and 
(ii) a tip with a slowly varying radius R = Rgmwit) but a 
fixed center of mass r^ip. Such a weak time dependence is 
necessary if one wants to approach the fluctuating brush 
surface without perturbing it too much during the initial 
stage of penetration . In both cases, initial parameters 
{R, rjVp(O) and v^^^ for case (i), and rtip, i?grow(0) and 

''^grow fo'^ <3ase (ii)) are also adjusted in order to produce 
the same quasi-static protocol with respect to the pene- 
tration depth etip(i) = etip (cf. Fig. [3]). 

The force F^ip required to push the spherical indenter 
into the brush at a given distance Ctip can be found in 
the limit etip <S^ H ^ R, these quantities being defined 
in Fig.[3l Depending on the ratio between H^ /R and the 
penetration depth e, one can distinguish two regimes: (i) 
the strong penetration regime (s), which is discussed in 
Ref4^, and (ii) the weak penetration (w) regime, which 
has also been considered in Ref.— . The latter regime is 
equivalent to the situation where an infinitely hard sphere 
is pushed against a film of modulus Gb. We only consider 
the latter regime and obtains the force i^tip^ 
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Figure 4. (Color online) Effective shear modulus of the brush 
as a function of the reduced grafting density p* . Results are 
plotted for the two types of tips discussed in the text. The 
theoretical prediction Gh = 8.5 ■ lO^^Pg^e/a'', from Eq. pUjI . 
is indicated by a dashed line. 
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Numerically, for the two types of tips, the average force 
exerted on the indenter is monitored and plotted for dif- 

3/2 

ferent penetration depths etip. The dependence on e^[ 
of Eq. ([TT|) is then only observed for the lower grafting 
densities because the lateral displacements of the grafted 
chains are less constrained by their stretching. When pg 
increases, the compressibility of the surface decreases as 
the density inside the brush increases. The short range 
lateral perturbation of the brush in the vicinity of the 
inserted tip may thus be supplanted by additional long 
range deformation modes, which are not taken into ac- 
count in the previous theory. 

In Fig. 01 the effective shear moduli Gh of brushes with 
varying grafting densities is shown for the two types of 
tips. The data are compatible with a quadratic depen- 
dence of Gh on the grafting density, p*, predicted by 
Eq. (Uni). Using v = 1/pi = 1.27a^ and Eq. ^, one 
obtains Gh = 8.5 • lO^'^pl^e/a^. The prediction qualita- 
tively agrees with the simulation data. Possible reasons 
for the quantitative deviations are the rather short chain 
lengths and the finite compressibility of the polymer liq- 
uid in the simulations, as well as the strong stretching 
assumption invoked in the analytical calculation. 



III. RESULTS 

When a liquid is in contact with a solid surface, the bal- 
ance of surface and interface tensions dictates the equi- 
librium properties of the liquid. Balancing the tensions 



Figure 5. (Color onhne) Density contour plot of a polymer 
drop (A*' = 10) wetting a brush (A'' = 40) for various compat- 
ibility parameters ebd ~ 0.3e (black) to tbd ~ 0.8e (orange) 
between the brush and the droplet. The reduced grafting den- 
sity of the polymer brush is pg — 13.1 and the droplet contains 
M — 3832 polymers. The dashed line represents the average 
free surface position of the unperturbed brush. The lifting-up 
of the contact line and the formation of a ridge is visible. In- 
set: Sketch of the procedure used to extract the static wetting 
contact angle, 6y, and the strength of the lifting-up, hr, here 
for a compatibility ebd = 0.6e. 



parallel to the surface at the three-phase contact between 
the solid, the liquid and its vapor, one obtains Young's 
equationi^ 7lv cos By + 7ls ~ 7VS: where 9y is the equi- 
librium contact angle, and 7lv, 7lSi 7vs a-rc the liquid- 
vapor, liquid-surface and vapor-surface surface tensions, 
respectively. If the surface is soft, then the forces that 
act at the three-phase contact line will also deform the 
surface and lift up the three-phase contact linei^"— . 

In our simulations, making the molecules of the brush 
and the polymer liquid slightly incompatible, we can in- 
dependently control the contact angle and the deforma- 
bility of the polymer brush. The brush is the softer the 
higher the grafting density Pg is. Therefore, quite low 
grafting densities are used in the following. The different 
wetting properties are explored as a function of ebd: a 
value Cbd/e = 1 leads to complete wetting, i.e. a thick 
polymer film spreads on the brush, while smaller values 
of ebd give rise to polymer droplets with a finite contact 
angle. We will refer to ebd as the compatibility parame- 
ter. 



A. Wetting of a cylindrical droplet 

Both, in experiment a^^'^^ and in computer simula- 
tions^"—, the determination of the apparent contact an- 
gle remains a delicate task. Commonly, for flat surfaces, 
an apparent, local contact angle is graphically extracted 
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Figure 6. (Color online) Equilibrium contact angle, Oy, be- 
tween the brush and the droplet as a function of compatibility, 
Ebd/e- Main panel: Contact angle for M — 1916 polymers of 
length A*' = 10 and the three lowest brush grafting densities 
p* = PgRc- The contact angle 9y slightly increases with p|. 
Inset: Contact angle for a fixed Pg = 9.5 and increasing the 
number of polymer chains in the droplet. The contact angle, 
9y, slightly decreases with the increase of p* at fixed tbd/e- 



at the contact line, from the value of the angle between 
an extrapolated tangent plane to the droplet's surface 
and the plane of the flat surface. Since we deal with de- 
forniable surfaces, the task is even more delicate, and 
we have opted for a carefully estimate from the den- 
sity contour plots, such as ones depicted in Fig. [5] The 
method of determining the contact angle is illustrated 
in the inset of this figure. In order to obtain an accu- 
rate estimate of the contact angle, the static properties 
of the brush/droplet systems were evaluated by collect- 
ing 10^ snapshots over simulation runs of 500r. In Fig. [SI 
the density contours arc plotted for a liquid droplet of 
M = 3832 polymers with A^ = 10 beads per chain. The 
drop is in contact with a brush surface of grafting density 
p* = 13.1, and contours are plotted for various compati- 
bilities 0.3 < Ebd/e < 0.8. 

As expected, it appears in Fig. [5] that the droplet de- 
forms its shape when the compatibility with the surface 
increases, leading to (i) a contact angle variation, and 
(ii) to the appearance of a lifting-up of the surface at the 
three-phase contact line. 

In Fig. [S] the dependence of the contact angle on the 
compatibility ebd is presented. Droplets with a fixed 
number of polymers and a degree of polymerization N — 
10 wet brushes at intermediate grafting densities. As one 
increases the compatibility, the contact angle decreases^. 
It also appears a slight dependence on pg, which becomes 
even smaller as one increases the grafting density (not 
shown), because at very high grafting densities the sur- 
face of the brush in a bad solvent is largely indepen- 
dent from the grafting density, pg. Similarly, also brush- 
droplet interface is rather independent of pg for dense 



Figure 7. (Color online) Dependence of the height h^ of the 
ridge at the three-phase contact line as a function of the poly- 
mer drop (A'^ = 10) and brush (A'^ = 40) compatibility ebd/e- 
The drops contain M = 1916 polymers, and low to inter- 
mediate brush grafting densities are considered. The dashed 
lines are fits using Eq. (|12|) with contact angles from the 
Fig. [6l Inset: wetting ridge profile for a M = 3832 poly- 
mer chains droplet with a chain length A'^ = 10. The droplet 
wets a brush of grafting density Pg = 9.5 with a compatibility 
Ebd/e ~ 0.7. The rescaled distance from the origin y* /a, is 
defined at the three-phase contact line. The Eqs. (|12p and 
(|13p are plotted, respectively with dashed and mixed-dashed 
lines, for 7d = 7lv, 7bv and Gb taken from the Tab. IT] Fig. [2] 
and Fig. |4l also respectively. The arrow marks the average 
brush height far away from the ridge. 



brushes. On the other hand, at low grafting densities, 
the decrease of the contact angle at fixed ebd reflects the 
intricate changes of the brush-melt interface in response 
to changing pg. We note that a value of Cbd/e ~ 0.6 is 
sufficient to achieve a contact angle of 90°, which is ex- 
pected to result in the largest deformation of the brush 
at the three-phase contact line. 

The deformation of the surface and the building of 
a wetting ridge at the three-phase contact line are dic- 
tated by a competition between the minimization of the 
droplet's surface and the deformation of the brush sur- 
face and brush-melt interface, respectively. In the inset of 
Fig. ini the reduced grafting density of the brush is fixed 
to Pg = 9.5, while the number, M, of polymer chains in 
the droplet increases. We observe that at fixed compati- 
bility, Ebd/e, the contact angle decreases when the size of 
the droplet increases. Hence, smaller drops seem to bal- 
ance surface, adhesion and elastic deformation energies 
in a different manner than larger drops. 

In Fig. [71 the height h,- of the lifting up at the three- 
phase contact line is presented as a function of the com- 
patibility Ebd/e. Data are extracted from the density con- 
tour plots, such as the ones drawn in Fig. \5\ Results are 
only shown for the A^ = 10 beads droplets that contain 
M = 1916 polymer chains, the trend is the same for 



the other droplet sizes and degrees of polymerization, as 
also noted in Ref J^. We consider small to intermediate 
brush grafting densities. Because the main deformation 
is localized at the edge of the brush-droplet interface, 
its strength is of the order of the monomer size a, and 
it slightly increases when the grafting density decreases 
and the surface becomes more deformable. We observe 
a weak dependence on the brush grafting density p* in 
particular because the deformation is rather localized in 
the vicinity of the three-phase contact line and the narrow 
brush-melt /brush- vapor interface. We have also observed 
(not shown) that the vertical displacement h^ at the con- 
tact line has a magnitude that depends on the droplet 
sizei^ii^ as ~ (7d/G'b)ln(i?d/cr), where 7d and Rd are 
respectively the droplet surface tension and its radius. 
This dependence has been checked for an estimate of R^ 
from density contours, as well as from a rough estimate 
for cylindrical droplets, i.e. Rd ex \/N x M. 

It appears that the height, h^, of the ridge increases 
when the compatibility increases to 0.8. This behavior is 
surprising because the contact angle decreases with Cbd/e, 
and the liquid- vapor interface tension acts more tangen- 
tially to the surface. In the limit Cbd/e ~> 1; the liquid- 
vapor interface makes zero (wetting) or a very shallow an- 
gle (autophobicity) with the brush surface. Starting from 
a droplet geometry with a large contact angle, thus, we 
expect a deformation of the contact line with an increase 
of the wetting ridge to a maximum value, followed by a 
decrease towards zero, as we increase the compatibility. 
Due to the softness of the surface that also slows down 
droplets spreading, we were not able to computationally 
access this time scale. 

The resulting wetting ridge of the soft surface has a 
height h^{y) that decreases from the three-phase contact 
line position. It has been theoretically proposed by Carre 
et ali^— and verified in experiments, that hy:{y) has an 
asymptotic behavior of the form: 



hr{y) ~ h,{0) In ( - ) , where h,{0) 



7d sin 0Y 
27rGb 



(12) 



where d represents the distance away from the three- 
phase contact at which no displacement occurs within 
the surface, 7d the droplet surface tension, and Gb the 
shear modulus of the deformable solid. The height at 
the origin, /ir(0), is taken as the reference position of 
the three-phase contact line. In the Fig. [3 we plot the 
this height as a function of the compatibility ebd/e, us- 
ing the associated contact angles from the Fig. [6l First, 
the theoretical prediction is found to underestimate the 
strength of the lifting-up of the contact line measured 
in the simulations. This prediction also depicts a maxi- 
mum at 9y = 7r/2, before decreasing at smaller angles; 
a trend that is not observed in our simulations, where 
/ir(0) continuously grows with the decrease of Oy- On 
the other hand, we find that the logarithmic dependence 
on the distance from the contact line is compatible with 
our simulations data, as shown in the inset of Fig. [71 



where the height of the ridge at the contact line is the 
one measured in our simulations. 

Eq. ([T^ predicts a logarithmic dependence that be- 
comes singular at the three-phase contact line. It has 
been recently proposed^ to model the shape in the vicin- 
ity of the apex by a cusp: 



h,{0)-h,{y) = ^^sme\y 
7bv 



(13) 



which then assumes that in the vicinity of the three- 
phase contact line, the surface elasticity does not affect 
the shape of the cusp. Such a behavior is plotted and 
compared to the ridge profile in the inset of the Fig. [71 
where a linear dependence is indeed compatible with the 
simulation at small distances from the apex, d < 5a. 



B. Wetting of a spherical droplet 

In this section, we study spherical, cap-shaped droplets 
in order to investigate the line tension on a soft, de- 
formable substrate. 



1. Line tension effect 

When one deals with spherical droplets, one also has 
to account for the curvature k of the three-phase contact 
line, which may have a significant contribution on the 
wetting properties of small enough drops. This effect is 
commonly represented via the line tension, r, that acts 
over the circumference length 27ra of the wetting sphere- 
cap shaped droplet, and modifies Young's equation ac- 
cording to: 



cos6'y = cos^oo K 

7d 



(14) 



with 7d the droplet surface tension with its vapor, and 9oo 
the contact angle of a droplet with zero curvature, typ- 
ically either a macroscopic droplet of a cylindrical one 
that has a straight contact line. The contact angles 9y 
are extracted from the equilibrium profiles of spherical 
droplets ranging from M = 984 to M = 30656 poly- 
mer chains of length A^ = 10. The drops wet different 
deformable surfaces of same grafting densities pg than 
the ones already used in the Sec. IIIIAI The brushes 



contain between M = 3969 , (L^ 



212cr) and 



M ~ 11881 , (Lx ^ Ly ~ 352(7) grafted chains of length 
N ~ 40, hence the lateral dimensions are large enough 
to avoid finite size effects. 

Results for cos Oy versus the curvature^ of the three- 
phase contact line k are plotted in Fig. [51 The dashed 
red line is the best fit that assumes a behavior according 
to the Eq. dHl), with 9^ = 78.61°±5.94° and the charac- 
teristic length scale r/7d = 18.96 ±6.040-. In the limit of 
large spherical droplets for which the contact line term is 
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Figure 8. (Color online) Cosine of the spherical droplet con- 
tact angle, 9y, as a function of curvature, k, at the three-phase 
contact line. The slope of the dashed red line corresponds to 
the characteristic length scale, r/7d ~ 19cr. The same kind of 
curves can be observed experimentallj"^ for nanometer-size 
sessile droplets. Inset: droplet radius _Ri vs the contact ra- 
dius K^^ , measured independently from the simulations. The 
blue line is a square-root fit to the spherical-cap approxima- 
tion, which is more accurate for large and slightly deformed 
droplets. This gives an estimate of the critical droplet radius. 
Re ~ 15.7(T, below which line tension effects become impor- 
tant. 



Figure 9. (Color online) B/L/V Interfaces in a minimalistic 
droplet/soft-surface model. The shapes are obtained by mini- 
mizing the Eq. (|15|l . in which surface tensions and brush elas- 
ticity are mapped onto simulation systems. An initial droplet 
of radius Rp < Re is considered, and the interface Hamilto- 
nian is minimized including or not the line tension term r/^d 
from the Eq. (|14|l . Results from the polymer model are also 
superimposed for the same size regime, contact angle and a 
brush grafting density p* = 13.1. Inset: a larger droplet 
of initial radius Rp > Re is considered, for a contact angle 
610 ~ 27r/3 and an effective grafting density Pg = 6.5. 



takes the form: 



not relevant, 9oo is then of same order than 6y measured 
for cylindrical droplets of equivalent surface contact ra- 
dius (inset of Fig. |6] for the largest system) . 

Therefore, one can relate the contact radius k~^ to the 
apparent droplet radius Ri measured independently. For 
large and slightly deformed droplets, the spherical cap 



approximation connects the contact radius 



to the 



square-root of the apparent droplet radius i?i. This is 
plotted in the inset of Fig. [S] in order to fit the measured 
data. Using the threshold value r/7d, one can estimate 
the critical droplet size. Re ~ 15.7a, below which the line 
tension may affect, at equilibrium, the wetting properties 
of spherical droplets^. In Fig. [SJ this corresponds to the 
shaded region. To reduce significant curvature effect at 
the contact line as well as dependence of the contact angle 
with the droplet size and brush/droplet compatibility, 
one then has to consider droplets with an higher radius. 
One way to quantify the line tension effect consists in 
using a schematic model that incorporates the free-energy 
cost stemming from the curvature of the contact line. To 
this end, we use an effective Hamiltonia n^°'^^ , which also 
allows us to qualitatively study the interplay between the 
surface and interfacial tensions, the line tension, and the 
softness of the substrate that dictates the droplet shape 
and the wetting ridge. The effective interface Hamilto- 
nian, H, accounts for the free energies of the droplet 's 
surfaces, the elastic properties of the deformable surface, 
and the line tension. For a three dimensional system, Ji 



n = Y,jiAi +tL + XVd + PHb 



(15) 



where the sum runs over the different interfaces / of 
area A/, L is the length of the three-phase contact line, 
and A is a Lagrange multiplier for the fixed volume con- 
straints of the droplet. The brush-like surface is mod- 
eled by vertical harmonic springs of stiffness k propor- 
tional to the grafting density, as already discussed in the 
Sec. IIIBI The stiffness is adjusted in order to approx- 
imate the one of our soft surfaces, and using Eq. ([TU| . 
we get k = 15Gb/(i?cPg)- This leads to the additional 
energy contribution Hb ~ J k/2{p — po)^dr in the to- 
tal Hamiltonian (|15p . which can be discretized in lo- 
cal densities Pi^^j, of spring segment lengths. Therefore, 

one can write Hb = J2i,^ii fcAT4oii/2 {pi^^^ - po) , where 
po is the initial spring segment lengths density in cells, 
AVcoii = (ij/coii ^ ^z^oii) ^'^"i t^*^ edge length L^^^i^ of the 
cell is half the spring equilibrium length. Other interfaces 
are also modeled by springs of stiffness inversely propor- 
tional to the spring equilibrium length. The Eq. (jlSp 
is then minimized, using a Monte-Carlo scheme, for dif- 
ferent initial droplet states 9o and radii Rp, and for a 
spreading coefficient S = 7bv — 7bl — 7lv that is inde- 
pendently tuned in the partial wetting regime, i.e. S < 0. 
In the main panel of Fig. [5J results from minimizations 
of Eq. P^ are plotted, in rescalcd units, for an initial 
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2. Wetting regimes 

In the following, we aim to evaluate how the charac- 
teristics of the droplet and the soft surface, can dictate 
the equilibrium shape of the droplet. We first estimate 
the free energy that can be gained, by a droplet, by de- 
forming the surface. We consider two limiting case, (i) 
wetting of a droplet over a hard surface of shear modulus 
Gb ^ oo, and (ii) wetting over a liquid surface (Gb -^ 0) 
and for which the surface deformation S is small com- 
pared to the droplet height h (see Fig. [10] for a sketch of 
the model). 

In the case (i) and using notations from the Fig. [10] 
with (5 — )■ 0, the balance of tangential forces at the three- 
phase contact line gives: 



AE, 



(') 



(7BL - 7BV) + 7LV ( 1 - -^ 



(16) 



Figure 10. (Color online) Two dimensional schematics of the 
spherical droplet model (L) at the interface between a sur- 
face (B) and the vapor (V) . Also drawn are the initial droplet 
radius Rp, its upper (-Ri) and lower (R) spherical cap approx- 
imation radii, the projected contact radius a, the penetration 
depth S and the visible droplet height h, as well as the two 
relevant contact angles 6 and 9' . 



where the 7ij are the surface tensions of the differ- 
ent interfaces drawn in the Fig. 1101 Using the spheri- 
cal cap approximation, the spherical cap radius writes 
Ri = (a^ + h^) /2h, and one obtains for the free energy 
cost: 



AEu 



(i) 



na 



7bl - 7bv + 7lv o , ,o 



(17) 



droplet of radius less than R^ and 6q = 7r/2, and for a 
value of 7bl that fulfills the Young's equation for 6*0- For 
quantitative comparison, rescaled data from MD simu- 
lation on the polymer model are superimposed, and for 
the same size regime, contact angle and surface softness 
of p* = 13.1 . For the minimization procedure, we have 
also considered two surfaces that mimic the ones we use in 
our simulations by adjusting the brush stiffness k accord- 
ingly. As in the Fig. [7] the strength of the ridge slightly 
depends on the surface compliance, at least in the range 
of stiffnesses we considered. Its strength and trend also 
roughly follow the ones given in the Fig. [7] for the current 
value of the contact angle, as well as compared to MD 
data on the polymer model. We also observe a decrease 
of the ridge when a line tension of same order than the 
one we measured in our simulations is included in the 
Eq. (jlSp . which therefore increases the equilibrium con- 
tact angle value, as also observed in Ref<^. In the inset 
of the Fig. in] the initial droplet radius Rp is larger than 
Rc^ while the initial angle 9q is equal to 27r/3, the surface 
softness mimics a grafting density p* = 6.5, and Eq. p^ 
is minimized with and without the hue tension term. In 
this droplets size regime, the line tension strength is sup- 
posed to not affect the equilibrium contact angle, which 
is indeed observed. 



In the case (ii), the lense-shaped droplet of the Fig.lTOl 
involves a free energy cost, when it wets a liquid, of the 



form AE, 



(vi) 



-na-^S 



lujh^ + 7bl^^ 



where the 



spreading parameter writes S = 7bv — 7bl — 7lv = 
7Lv(cos0y — 1). Therefore, in the spherical cap approxi- 
mation, the difference AE ~ AEi^a) — AEi^i) provides an 
estimate how much free energy can be gained by deform- 
ing the surface, and one obtains: 



AE = -K { 7bl<5 + 7lv^- 



3a2 



h^ 



(18) 



Then, we estimate the free energy cost of deforming a 
soft surface that is mimicked by a polymer brush, thus 
with the constraint of attachments on one chain ends. 
Using notations from the Fig. [TU] this one involves the 
formation part of the free energy of the drop, and the 
one of a round AFM tip of radius R that gives rise to the 

Eq. (dH), namely Uei = (32/15)Gb#^^<55/2. For the for- 
mation of the drop, the free energy is split in the energy 
cost of displacing an area tto^ of B /V interface, and the 
one that comes from the upper (i?i) and lower cap radii 
(i?) and their related interfaces. The final free energy is 
then written in the spherical cap approximation, which 
gives: 
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AF 



32 

15%/2 



G'b(a' 



^2Nl/2^2 



Tra 5- 



^ (7Lv/i' + TBL*?') , (19) 



where Gb is the shear modulus of the soft surface. Thus, 
Eq. p^ approximates the rubber surface as a hquid-hke 
surface that has an elastic response under deformation. 
It combines both the droplet deformation under wetting 
and its effect on deforming the surface. 

By comparing Eqs. ([T5|) and P^ . we note that the 
asymptotic limit AF/ AE ^ 1 corresponds to a hard sur- 
face for which the free energy cost is prohibitive, and for 
which the Young's equation dictates the droplet shape. 
The opposite limit AF/ AE — >■ corresponds to a very 
cheap cost of deformation, and the lens-shaped droplet 
wets a liquid. At intermediate values, AF/ AE < ±1, a 
concurrent coupling appears between the elastic proper- 
ties of the drops and the surface, and therefore, a more 
intricate behavior dictates the droplet shape. One also 
should note that for AE > it is more favorable for the 
drop to wet a solid, while for AE < 0, it costs more. 

In order to extract a cross-over length, and follow- 
ing Refs4^i^ for a small perturbation of the brush 
surface, we approximate the elastic free energy Uei ~ 
AbGb{S^a^)/{Hf = AbGya^/{Hf, where {H) is the av- 
erage brush height and A^ is given in the Eq. (j20b[) . Us- 
ing notations of the Fig. [TUl we also rewrite Eq. p^ in 
a more compact and completely equivalent fashion: 



AF = AbGbaV(^)' + AdiBY^a^ (20a) 

Ad = 7bv [S + 7LV tan^ {9/2) + jbl tan^ {9' /2)] 



and, Ab = (37r/2) tan^ {9' /2) 



(20b) 
(20c) 



As we are seeking for a preferable droplet scale for 
which wetting over a deformable surface is stable and 
favoring, we then minimize the free energy AF with the 
constraint that the droplet volume is constant. For this 
purpose, one defines the Lagrange function: 



L (a, 0, 9', A) = AF (a, 9, 9') - A [Va {9, 9') - V^] (21a) 
and, Vd [9, 9') = A^a^ (21b) 

= (7r/6){3 tan (9/2) + tan^ (61/2) 
+ 3 tan {9'/2) + tan^ (0' /2)}a^ . 



where, in Eq. (|21ap . A is a Lagrange multiplier and Vp = 
(AttR'^) /3 is the initial droplet volume, using notations 
from the Fig. [TUl The principle of least energy and the 
condition of dL/d9 = yield to 



A = 4a'^7Lv tan (61/2) cos^ (61/2) 



The condition dL/da = 0, in which the result from 
Eq. (22) is inserted, gives: 



= tan (61/2) cos^ (61/2) . (23) 



AbGho'^ 7rArf7Bv 



2A„7Lv(^)3 67Lv^. 



Finally substituting Eq. p2|) in the condition dL/d9' = 
0. one obtains: 



tan (6172) cos^ (6172) 



3GbQ'^ 7BL 

27Lv(-ff^> 7LV 

X cos^ (61/2) 



tan (61/2) 
(24) 



In view of Eq. (P^ , one can note that the right-hand side 
is of 0(1), which can be fulfilled on the left-hand side by 
requiring the condition: 



AbGh (a*) TrAd-jBY 



2(if3) 



6 



(25) 



where a* appears to be a characteristic scale. Using the 
relation (H) ^ vNpg = vp*N/Rl discussed in Sec. IIIB[ 
one finally arrives to the expression of a cross-over length 
scale at which the free energy cost of deforming a soft 
surface with respect to the droplet shape is minimal: 



--^1) 



3/4 



7i"Ad7Bv {pIY 
3AbGh 



1/4 



(26) 



(22) 



where A^ and R^ arc respectively the number of 
monomers in a brush chain, and the end-to-end aver- 
age radius for this length A^. This scale marks the bor- 
der between a domain where droplets of radius lower 
than ^ have an equilibrium shape dictated by a surface 
tension/elasticity coupling, to another domain for larger 
droplets where elasticity dominates. 

To evaluate Eqs. (HH)), (HH), and 1^ over different 
regimes, the size of the droplets was controlled by increas- 
ing the number M of polymer chains inside the droplet, 
i.e. Rp ex M^/"^, while different wetting regime were con- 
sidered by either tuning the brush/droplet compatibility 
fbd/e and the brush grafting density p*. In any case, only 
droplets of radius Rp > Re were simulated. 

In Fig.[TTl we have plotted the ratio AF/ AE as a func- 
tion of ^/a, for data obtained from simulations. Such a 
ratio gives a complete map of trends for the wetting of 
droplets over either a hard, liquid, or soft surface, respec- 
tively described by AFj,), AFjj.j), and AF. Note that 
the soft and liquid surfaces, within this model, mainly 
differ by the elastic energy term of deforming the sur- 
face. The different wetting regimes are labeled as follow, 
A: for the droplet, wetting a hard surface is better than 
a soft surface, while the equilibrium shape is driven by a 
surface tension/elasticity coupling, B: the same wetting 
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energy. 



IV. CONCLUDING REMARKS 



' [^a](p*g'M) 



Figure 11. (Color online) Cross-over diagram between the 
free energy cost AE of deforming a solid surface and a liquid 
along the lines of Eq. (|18|l . and the free energy cost AF of 
deforming an elastic soft surface from Eq. (|19|l . experienced 
by a spherical droplet. The ratio AF/ AE is plotted as a 
function of the cross-over ratio ^/a from the Eq. (|26p . which 
define four different wetting regimes {A,B,C,D} as discussed 
in the text. Inset: scaling plot which points out the relevant 
role of the adhesion energy Wa for droplets larger than ^, 
which experience a coupling between surface tension effects 
and brush elasticity. 



regime is favorable but surface tension dictates the equi- 
librium droplet shape, C: idem for the shape but it is 
more favorable for the droplet to wet a soft surface, and 
D: the same wetting regime might be preferred whereas 
the droplet shape is achieved in a coupling between sur- 
face tension/elasticity. 

First we note that droplets smaller than S, tend to be 
located in regime B for larger brush grafting densities 
and in C otherwise. The shape of the smallest drops does 
not strongly differ between a hard or a soft substrate, the 
minimization of the surface tension dictates the behav- 
ior. A balance is achieved between surface and adhesion 
energies, leading to the Young-Dupre equation^. Sec- 
ond, one observes for droplets of size larger than ^ that 
increasing p* mainly drives the systems from regime D to 
A: the higher p* is, the larger must be the compatibility 
Ebd/e in order to for the melt to penetrate into the brush 
and reach a more favorable "soft-wetting" state, which 
is mainly driven by the competition between the brush 
elasticity term in Eq. (|19p and the brush/ liquid energy 
strength in the Eq. (fT8|) . Because one can expect the 
adhesion energy Wa to be proportional to Cbd^'^, one 
can confirm this trend by rescaling the ratio AF/AE by 
7d/Wa, where the droplet surface tension 7d = 7lv- Re- 
sults are summarized in the inset of Fig.[l] where, despite 
of the statistical error of the simulation data, the data are 
compatible with a scaling behavior. Thus, the behavior 
in regimes D and A is mainly governed by the adhesion 



Molecular dynamics simulations were performed in or- 
der to study the properties of the solid/liquid/vapor in- 
terfaces between viscous droplets and a deformable sur- 
face. First, particular attention has been paid to char- 
acterize the soft polymer brush surface. To this end, 
an implementation of an AFM-like numerical experiment 
has been developed, which has allowed us to extract the 
shear modulus of the brush as a function of the graft- 
ing density for the brush. The simulation data agree the 
ones theoretically predicted by Fredrickson et al^ and 
the experimental work of Fujii et al^. 

The wetting properties of the surface are determined 
by its softness and the compatibility between the droplets 
and the surface. The former can be controlled by varying 
the grafting density of the brush, the latter is strongly 
influenced by the compatibility between brush and liq- 
uid. Density contours have shown a deformation of the 
wetting ridge near to the three-phase contact line. The 
strength of this lifting-up was found to increase with the 
compatibility between brush and liquid, as well as with 
the capacity of deformation of the surface. 

Finally, the wetting properties of spherical droplets 
that are in contact with the same kind of deformable 
surface have been also investigated. Different drop sizes 
and brush grafting densities were considered. We esti- 
mate the line tension and the concomitant critical droplet 
radius, above which line-tension effects may not be deci- 
sive. Then, minimizing an effective interface Hamiltonian 
with effective parameters that mimic the thermodynamic 
properties of the simulated systems, it has been possible 
to capture the lifting-up of the contact line, and to re- 
cover the line tension effect for droplets of sizes larger or 
smaller than the previously mentioned critical radius. 

The interplay between the contact properties of the 
spherical drops and the deformable surface and its im- 
pact on the droplets shape deformation has then been ad- 
dressed within a similar effective interface Hamiltonian. 
As the surface/droplet compliance can be tuned inde- 
pendently to the surface softness, and also because one 
can deal with different droplet sizes, it has been possible 
to cover a wide range of wetting and deformation behav- 
iors. Two distinct asymptotic regimes have emerged. For 
droplets smaller than a critical size, the contact proper- 
ties and its impact on the drop shapes are mainly dic- 
tated by the balance between the adhesion and surface 
energies. In this limit, the drops equilibrium properties 
are well characterized by the Young-Dupre formalism. In 
an other limit, for larger radii, the balance of surface ten- 
sion and brush elasticity dictates the wetting behavior, 
and adhesion energies finally determine the equilibrium 
droplet shapes. 

The same kind of size effect has been shown to have 
a significant role on the dynamical properties of rolling 
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droplets over hard surfaces^. It is then of interest to 
extend the current study to the steady state motion of 
drops over a soft surface, first to evaluate the contribution 
of the hfting-up effect to the dissipation of the moving 
contact line, and then to address droplet-size effects of 
the dissipation mechanisms. Further investigations are 
also required to fully characterize the slippage of the soft 
surfaces when complex fluids low past it. 
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